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ABSTRACT 



Far field computational boundary conditions for 2D internal flow problems 
are developed from analytic solutions of the linearized Euler equations. The 
Euler equations are linearized about a constant pressure, rectilinear flow 
which nay have streamwi se-normal variations in temperature and velocity as a 
result of entropy production in the nonlinear computational region. The bound- 
ary procedure can be used with any numerical Euler solution method and allows 
computational boundaries to be placed much closer to the nonlinear 
region of interest. 
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I. INTRODUCTION 



Numerical solution procedures for nonlinear fluid dynamic equations 
usually use one or more artificial computational boundaries located at some 
distance from the primary region of interest in order to limit the physical 
domain to finite size. If the flow crossing such a boundary (either inflow or 
outflow) is subsonic, then some type of computational boundary conditions must 
be imposed which simulate the influence of the true far field conditions at 
infinity. These boundary conditions must be such that waves crossing the 
boundary do not produce erroneous reflections back into the computational 
field to degrade the calculations. It is generally acknowledged that simply 
imposing free stream conditions (or conditions at infinity) at computational 
boundaries is usually inappropriate. Standard practice has consisted of 
locating the boundaries quite far from the region of interest in an attempt to 
simplify the boundary condition models and minimize any effects of 
inconsistent modeling. The net effect is a significant increase in the number 
of grid points required for an accurate flowfield calculation. 

A boundary modeling procedure for two-dimensional internal flows is 
presented which alleviates the difficulties mentioned above and also allows 
the computational boundary to be located much closer to the nonlinear region 
of interest. The procedure is limited to steady, inviscid flow, although the 
flow can be rotational. It represents a logical extension of the so-called 
characteristic boundary conditions commonly used with inviscid numerical 
solution methods. Extension to axisymmetric or three-dimensional flows is 
straightforward. 

The analysis presented here is based on the Riemann variable formulation 
of the Euler equations given in Reference 1. This represents a natural 
starting point because the characteristic (or zero-order) boundary conditions 
mentioned above are expressed in terms of Riemann variables. The equations 
are linearized about a constant pressure, rectilinear flow condition, which 
truly represents conditions at infinity. These linearized equations are 
assumed applicable in the far field region beyond a computational boundary. 
Within the nonlinear computational domain, strong entropy-producing (i.e., 
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rotational) effects can exist which create variations in density, velocity, 
etc. in the far field in the streamline-normal direction which are not 
necessarily small perturbations. Such variations are modeled in the present 
analysis . 

The linearized equations are solved analytically using Fourier analysis 
techniques as outlined in Reference 2. These solutions are coupled to the 
nonlinear numerical solution to provide a smooth transition across the 
boundary to the true far field conditions at infinity. The coupling is 
accomplished by the boundary conditions. The underlying principle is that the 
streamwise variations of both upstream and downstream running waves should 
decay to zero at infinity. These first-order boundary conditions provide 
distributions of flow quantities to be imposed along the boundary, not 
constant conditions. They represent a logical extension of the zero-order (on 
characteristic) boundary conditions. Furthermore, the boundary analysis can 
be coupled with any inviscid numerical solution method. 

The boundary condition analysis has been applied to two-dimensional 
duct flow and to cascade flow, where conditions are periodic in the direction 
of the blade row. For duct flow, both isentropic and non-isentropic boundary 
conditions are derived. Only isentropic results are given for cascades. 
Extension to non-isentropic cascade flow can be carried out by following the 
procedure used for duct flow. 

Numerical results are presented for both isentropic and non-isentropic 
duct flow. Results obtained using the first-order boundary conditions are 
compared with those using the zero-order boundary conditions. Numerical 
results for cascade flow will be presented in Reference 3. It was found thas 
the size of the computational field and associated number of grid points 
needed for the nonlinear numerical solution could be reduced significantly b 
using the new first-order boundary condition procedure with no loss in 
numerical solution accuracy. The reduction in number of grid points was as 
much as 50 percent in some cases. The additional computational effort 
required by the new boundary procedure is small (less than 10 percent) so th 
a significant saving in overall computational effort was realized. A large 
portion of the gain is due to the sizeable reduction in the physical extent 
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the computational field and the fact that fewer solution iterations are 
required for information to propagate between the upstream and downstream 
computational boundaries leading to more rapid solution convergence. 
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11. PERTURBATION EQUATION DERIVATION 



The system of two-dimensional, steady, linearized Euler equations which 
describe first-order spatial perturbations from a constant pressure state will 
be derived in this section. A Riemann variable formulation taken from 
Reference 1 will be used because of its close relationship with the character- 
istic (or zero-order) boundary conditions commonly used in numerical solution 
of the nonlinear Euler equations. 

The two-dimensional form of the Euler equations is (Reference 1) 



di 



aR . .s aR _ . v-l . , r 
St + <« ' a) SI - + 2 a (s 



Y-l' 



ae + ae = _ a£ aP 

at q as yq an 



as 

at 




o 



(1) 

( 2 ) 

(3) 

(4) 



Velocity magnitude and speed of sound are denoted by q and a, respectively, 
and P is the logarithm of pressure. The Extended Riemann Variables Q and R 
are defined as 

Q = q + aS (5 
R = q - aS 



while the modified entropy is defined in terms of pressure p and density p as 

S = [2y ” log (p/p Y )] ( 6 ) 

The flow angle is 0, time is denoted by t, and local distances along and 
normal to the streamline direction are denoted by s and n, respectively. 
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For steady flow the analysis can be greatly simplified by defining a new 
dependent variable 



The local Mach number is denoted by M. Equation (8) is obtained by 
subtracting equations (1) and (2). Equation (9) is obtained by adding 
equations (1) and (2) and integrating. The constant of integration, which is 
proportional to stagnation temperature, can be set to unity by proper choice 
of non-dimensional izing quantities. The simplified form of the steady Euler 
equations is then 



According to equation (12) entropy remains constant along streamlines. 

In regions of the flowfield where nonlinear effects are weak, the flow can be 
treated as a perturbation to a constant pressure, rectilinear flow. Such 
regions occur near and beyond far field computational boundaries. The depend- 
ent variables in equations ( 10 ) - ( 12 ) can then be expanded in asymptotic series 



T = Q - R 



(7) 



and replacing equations (1) and (2) by 



(M2 - 1) U + (Y-l) 9 H S |i = 0 

a 2 + 111 q2 = i 



( 8 ) 



(9) 






( 10 ) 




( 11 ) 




( 12 ) 



a 2 + III q 2 = i 



(13) 



T = T + T, + T, + . . . 
» 1 Z 



s = s + s. + s, + . . . 

CO 1 2 



( 14 ) 



0=e + e. + 0 O + . . . 

CO 1 2 
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The flow direction at infinity is assumed constant and denoted by 0 ; the 

oo 

perturbation quantities T. , and 0 ^ vanish at infinity. Entropy variation 
is not excluded so that the flow can be rotational. Furthermore, entropy 
variations can be strong (i.e., not small perturbations) so that S and T are 

oo oo 

not necessarily constant, but may vary normal to the streamline direction. 

Note that T depends only on S (because p is constant). 

oo oo oo 

Consistent with (14), spatial derivatives in equations (10)-(12) can be 
approximated by 



f^r = cos ©„ (l~r 0 i)|^ + cos 6 oo ( T+ ® 1 ) 3 ^ + ... 



(15) 



= -cos 0oo (t+0l)|- + COS 0„ (1-T0l)|y + .. 



where x and y are reference Cartesian coordinates, 0 is measured from the x 
axis, and 



t b tan 9 



(16) 



If expansions (14) and (15) are introduced into equations ( 10) - ( 12 ) , the 
resulting first-order perturbed Euler equations are 



, m 2 _ nr £ll . Hi 

^ ^ax T ay 



SToo 8T* 

9 l •' 3F ‘ 5T )] * 



( Y - i ) q« Moo Soo [ 



801 

ay~ 



801 

' 0^ = ° 



(17) 



Hi M 2 C r^l 4. ■ 1 r aT l __i . r - ~ 

2 Moo s„ [g x T ay J 2a„ * 8 y T 8 x T 0 



8 T 1 Ti 8T 0 

ax TI *-8y" 



aToo n , 

' 5T ]) + 



r± {s . r!il - , Hi] + r!!: 

2 l 3y ax j So,, l 8 y 



8 Soo 



(18) 



8Si 8Si 8 So, 

8^ + T 8T + 0 1 [ i7 



a Soo, 

T a = 0 



(19) 
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Velocity, speed of sound, and Mach number at infinity are denoted by q , a 

oo oo 

and M , respectively, and may vary normal to the streamline direction. The 

oo 

fact that all dependent quantities are convected unchanged along rectilinear 
streamlines at infinity has been incorporated into the above equations. 

Asymptotic expansions of the Riemann variables Q and R can also be 
defined as 



Q = Q„ + Qi + Q 2 + • • • 

R = R + R, + R_ + ... 

00 1 z 



(20) 



Using the definition (7) and the expansion (14) for T, it follows that 



T 

00 



R 



00 



T 



1 





( 21 ) 



Introducing the expansions (14) and (20) into the algebraic equation (13) 
gives the first-order relationship 

M* (Q 1 -R 1 ) + ^ m£ (Q 1 +R 1 ) - 2q„ Si = 0 (22) 

This will be used later in Section IV where the boundary conditions are 
derived. 
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III. SOLUTION OF FIRST-ORDER EQUATIONS 



Solutions of the first-order equations ( 17 ) - ( 19 ) are developed in this 
section for two-dimensional duct flows confined by parallel walls and for 
infinite cascade flows. Duct flow solutions are obtained for both isentropic 
and non-isentropic conditions. The cascade flow solution is provided for 
isentropic conditions. 



Duct Flow 



For two-dimensional duct flow between parallel walls, the far field flow 
angle 6 (and consequently t) can be set to zero by aligning the coordinate x 

co 

with the duct axis (see Figure 1). The width of the duct can be taken as 
unity without loss of generality. Equations ( 17 ) — ( 19 ) then reduce to 



? 3Ti 3T«, 30i 

(Mi - 1 ) [g~ + 01 gy] + (y-D q~ K, s„ — = 0 



ay 



(23) 



301 i 3Ti Tj 3T 0 



Si 3 S 0 



2 Qco Mo ° S “ 3x A + 2 [ 3y A " T* 3y ^ + a °° S„ 3y 



+ ( Soo - -§y) 



2 „ 9S 1 



(24) 



ay 



= o 



3Sj 

JT 



+ 01 




(25) 



Isentropic Conditions - According to the scaled entropy definition (6), the 

2 

far field entropy S is -=y. Equations (23)-(25) then simplify for isentropic 

oo 

flow to 



p 2 8Ti _ aej. _ 

2q 0 cM 00 3x 3y 



(26) 



30j i 3Tj 
3x 2q a> M 00 3y 



(27) 
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where 



p = V 1-M 2 

oo 

Note that q and M are constant. 

oo oo 



( 28 ) 



Equations (26) and (27) can be solved by separation of variables by 
assuming 



T l = 2 i. «. F(y) 

0! = e Ax H(y) 



(29) 



where X is the unknown separation constant. Equations (26) and (27) then 
reduce to the ordinary differential equation system 



p 2 A F - H' = 0 
AH + F' =0 



(30) 



with boundary conditions 

H ( 0 ) = H( 1 ) = 0 (31) 

Primes denote differentiation with respect to y. Eliminating F gives the 
single equation 

h" + A 2 p 2 H = 0 (32) 

which has the solution 

H = kj sin Apy + cos Apy (33) 

To satisfy the boundary conditions (31), the constant k^ must be zero and 

A = + “ (n = 1, 2, ...) (34) 
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Using these results, the general solution of the system (26) and (27) can 
be written as 



Tl 

01 



= z Km 

l 



2qooMoo 

P 



cos nuy 



sin niry 



nirx/ p 

e + 



(35) 



Z K 2n 
1 



-2q 00 M 0 



cos nny 



sin nny 



-nux/p 



where Kj n and K^ n are arbitrary coefficients. 



Non-Isentropic Conditions - Although the system ( 23 ) - ( 25 ) is linear, it has 
non-constant coefficients since is independent and varies with the 
streamline-normal direction y. The thermodynamic variables are related by 



Hi Hi (_1_ . s ) 

2y 2 S-l ic ° ; 

a«> = Po e 



where p is the constant pressure at infinity. 



(36) 



To achieve an analytic solution, a second level of linearization can be 
introduced by defining a new variable o(y) by 



S i -^r (l-o) 

oo Y”1 



Since o is generally small, equation (36) can be approximated by 



(37) 



^00 



Hi 

P„ 2y (l+o) 



(38) 
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An approximation for is then provided by equation (9). Reference 
quantities defined by 



Hi 

S. = p~ 2y 

q« = ^5 
= 



( 39 ) 



can be introduced to simplify the notation. The approximate thermodynamic 
relations then become 



a® = (1 + o) 

q. = a- (1 - —) 

Y 1 M 2 

oo 

(1 - -It =§) 

4 - 

Too ~ a oo 

Y"1 

All quadratic and higher terms in a have been neglected. 



(40) 



Introducing the above approximations (40) into the governing equations 
( 23 ) - ( 25 ) gives 



/ P 2 . _4 o> a ^l n / . 4 o ' a9 l _ n 

y- 1 a2 ax " 2qooM ~ (1 y- 1 q2 J ay “ 0 



(41) 



On M M JL _£\ £^i 4. !li - o; !_/ C \ 

2qooMo ° (1 ■ Y-1 q2 ax ay 2a ~ 3y (oSl) 



(42) 



!li - _2_ do 
3x y - 1 dy 1 



( 43 ) 
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The parameter p is defined by equation (28) with M substituted for M . 

oo oo 

Equation (43) suggests that is of higher order than or 0^. If so, then 
the RHS of equation (42) can be set to zero and equation (43) becomes 
decoupled. Validity of this assumption can be verified later by evaluating 
numerical results obtained from this analysis. The final system of governing 
equations is therefore 



(P‘ 



0Ti 



00! 



H Ov u ' X fj. i . 4 O, UV 1 

y-1 :2* ax " 2qooMo ° (1 ' y-i o2 ay 



(44) 



- - 4 o 301 0Ti 

ZqJLU-^jf) aF + aT = 0 < 45 > 

'oo 

After these equations are solved for Tj and 0^, the perturbation variable Sj 
can be obtained from equation (43) by quadrature. 

As for the isentropic analysis presented above, a separation of variables 
solution approach is again appropriate. Analagous to equations (29) assume 



T, = 2q m M e Ax F(y) 

0j = e Xx H(y) 

Equations (44) and (45) then reduce to the ordinary differential equation 
system 



H ' = [XF + sl H ' 5 



M 2 



F' + XH = -Ar ^ AH 
y-1 a 2 



with boundary conditions 



(47) 



H ( 0 ) = H( 1 ) = 0 



(48) 
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Eliminating F from equations (47) gives the single equation 



H" + A 2 p 2 H = -§=■ ^ Co(H" + A 2 p 2 H + AM 2 F') + o'(H' + Am£ F)] (49) 

oo 

This equation is in a form suitable for solution by iteration. The lowest 
order solution (i.e., first approximation) satisfies 



II 

H 



+ A 



p 2 H = 0 



F 




I 



(50) 



Using these results in the RHS (denoted by <|i) of equation (49) produces 

H + A 2 p 2 H = — ^-r • -l- [o H + m£ o H ] = <J> (51) 

P*-q*- 

r oo 

The solution of this equation provides an improved second approximation for 
the solution of equation (49). Repetition of this iterative process generates 
ever-improving approximations. Only the solution for the second approximation 
will be given here. Whether or not it is sufficiently accurate can be 
assessed by evaluating the numerical results presented later. 

The solution for the first-approximation equations (50) was obtained in 
the above isentropic analysis as 

00 

H = E A n sin nny ; A = + (52) 

1 " P 

where the A are arbitrary coefficients. Note that the iterative solution 
procedure used here is consistent in that the isentropic solution is recovered 
for vanishingly small o. Using the approximation (52) the RHS of equation 
(51) can be written 

4> = — r [o' l n A n cos nny - n MT, o l n 2 A n sin nny] (53) 

Y_1 p 2 q£ 1 1 
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Equation (51) can be separated into component equations 



" 2 2 
H + A p H = 



n n r n n 



(n = 1, 2,...) 



(54 



where <t> n represents the nth component of the source term <f> defined as 



4 unAfl i -~ 

4> n = “f ~I? [o cos nuy - u Hf o n sin nuy] 



^ p 2 a ! 



(55 



For consistency, the separation constant A n must also be expanded as 



A n =±(^+6A n+ •••) 



(56 



so that equation (54) becomes 



2 2 



H + n u H =4) 2upn 6A A sin nuy 
n n n n n 



(57 



The solution of this equation which satisfies the boundary conditions (48) i 



1 r y 

H n = K n sin nuy + p 6A n A n y cos nuy + _i J <|> n sin nu(y-q) d n 



(58 



where K n is an arbitrary coefficient and 



Si " = srs; S ♦" sin dn 



(59 



This solution is valid regardless of the choice of sign in equation (56). 



The solution for the nth component of F can be obtained from the system 
(47) as 



K n ,, 4 o 



F n * + [t (! ' 3TT cos niT y " 6x n A n V sin niT V 

p t -i 



( 6 ( 



i f y 

+ — T \ ♦n cos ni, (vn) dn] 
nup i - 
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The choice of sign corresponds to that of equation (56). 



Using the above results, the approximate solution of the system (44) and 
(45) can be written as 



Tl 

6l 



I K ln 
1 



2c ^~ / , 4 o \ 

“ (1 ' ~1 tfS) C0S 



♦ l K 2n 
1 



T 

sin nny 






e^ + “n)x 



( ^ Ox 

(1 - — Y ~ 2 r 2 ) cos nny 



P q 



sin nuy 



oo 


-2qJi f n 


/ DlT , r _ , ® 

e ( a 6A " )x + z 


2qoAo f n 


l 


P 


P 


1 


h n _ 


1 


^n _ 



/ n Tf . r . \ 
e -(— + 6A n )x 



e'^ + 5A "> X 



(61) 



where 



1 r y 

f n - 3 <5* n A n y sin miy - — V <J» n cos nn(y-ri) dn 

J o 

1 r y 

h n s p 5A n A n y cos nny + — J <J> n sin mr(y-n) dq 



(62) 

(63) 



This solution reduces to the isentropic result (35) when o is zero. Using 
this solution, equation (43) can be integrated to obtain the entropy 
perturbation S^. 



Cascade Flow 

For isentropic cascade flow, equations (17)-(19) reduce to 

3 2 r 3T l . 3T 1-, "1 . "1 _ n 

+ 1 - w + ’ sr - 0 



(64) 
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90 1 30 1 . 1 r 0T l 9T 1-, _ n 

+ T w 2^1 r c ^r ' T - 0 

oo oo J 



ax 



(65) 



If the (x,y) coordinates are chosen such that the flow is periodic in the y 
direction (see Figure 2), then a separation of variables defined as 



Tl = 29 A F(5) e 



lrniy 



0j = H(0 e 



muy 



( 66 ) 



£ 5 innx 

can be used. The blade spacing has been taken as unity. Equations (64) and 
(65) are thereby transformed to the ordinary differential equation system 



3 2 (F' + tF) + tH' - H = 0 



H' + tH - tF' + F = 0 



Primes denote differentiation with respect to £;. 



The solution of the system (67) has the form 



(67) 



F = ae 






H = e A ^ 



( 68 ) 



where the eigenvalues are 



A = 



t ± i p ( 1 + t^) 

oo 

p2TT2 



and 



(69)1 



A + T __ 7 1 

a - r - + — 

Ax - 1 p 



(70) 



Using these results, the general solution of the system (64) and (65) can be 
written as 



16 



Tl 

91 



where 



I K ln 



E K 2n 



-2qJ^J * 

P 

1 

2qooH»i " 
~ 

1 



irnifAjx 



inir( A 2 x 



+ y) 
+ y) 



+ 



x i = 

x 2 = 
Note that 



t + ip(l + T^) 

00 

T - ip(l + t2) 

p2 + x 2 

Aj and A^ are complex conjugates. 



(71) 



( 72 ) 
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IV. BOUNDARY CONDITION DEVELOPMENT 



Examination of equations ( 1 ) - ( 4 ) shows that at a subsonic far field 
computational boundary there are three downstream-running waves and one 
upstream running wave. Therefore, the information available from the 
numerical solution is not complete and differs at upstream and downstream 
boundaries. The information lacking is provided by the boundary conditions. 

If the flow is supersonic, all waves are downstream-running and specification 
of boundary conditions is straightforward. 

Far field computational boundary conditions (subsonic) are developed in 
this section based on the linearized Euler solutions obtained in the previous 
section. These solutions are assumed valid in the region beyond the computa- 
tional boundaries where nonlinear effects are small. Within the computational 
boundaries the full nonlinear Euler equations must be solved numerically. The 
boundary conditions provide for a smooth coupling of the nonlinear and linear 
solutions so that the true conditions at infinity can be imposed. 

The three specific cases analyzed in the previous section will be 
addressed in this section. 



Duct Flow 

Both upstream and downstream duct flow boundary conditions will be 
derived for isentropic conditions. The upstream boundary analysis is valid 
even if there is nonlinear entropy production downstream of the boundary 
within the computational domain. For non-isentropic conditions only 
downstream boundary conditions will be derived, since non-isentropic 
downstream flow is the more common situation. 

For isentropic, constant pressure flow in the upstream or downstream 
regions, the far field Mach number M^ is determined implicitly from the mass 
flow per unit area w by the relation 
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I 



1 1+Y 

M.U + ^ Mf] 2 1 ''' = w (73) 

The mass flow is usually known (or can be calculated) for a given duct flow 
problem. The associated speed of sound and velocity are 



cl oo 



ill 

M 1+ Y 

(—) 

w 



9oO ^OO Moo 



(74) 



Therefore, the far field quantities and R m appearing in expansions (20) are 



Qoo = q« + 



2 

Y" 1 



a« 




(75) 



p 

The far field entropy has been set to Equations (21) relate Q and R to 

the variable T. 



For non-isentropic, constant pressure flow in the downstream region, 
equations (37), (39), and (40) give the approximate far field relationships 



Qoo doo 




2 a oo 




o 



D ; 2 - 2 a °° 

Koo ~ doo " V _1 a oo " i - 0 

Y 1 Y 1 Moo 



T 4 - 

loo - a 0 



Quadratic and higher terms in o have been neglected. 



(76) 



Isentropic Conditions - At a computational boundary (assumed located at x=0) 
the perturbation flow variables can be represented by the Fourier series 
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(77) 



oo 



0j = I A n sin nuy 
1 



co 



Qj = l B n cos nuy 
1 



(78) 



oo 



Rj = I C n cos niry 
1 



(79) 



oo 



oo 



Tj = l E n cos nuy = l (B n -C n ) cos nuy 



(80) 



1 



1 



The boundary conditions will be developed from relationships between the 
Fourier coefficients A^, B n> and C^. The absence of modes corresponding to 
n=0 for Qj, R^ and T^ is related to the fact that these first-order 
perturbations must vanish at infinity. Further discussion of this topic is 
presented below in conjunction with the boundary condition development for 
cascade flows. 

The general solution for linearized isentropic flow is given by equation 
(35). Applying this solution at x=0 and using the series expansions ( 77 ) - ( 80 ) 
it follows that 



For the region upstream of the computational boundary (i.e., x<0), the 
exponentially growing component of the solution (35) can be suppressed by 
forcing K^ n to be zero. This requires the Fourier coefficients to be related 




(81) 




Therefore, 
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An = 2qJ^ (B " ' C " } 



(84) 



The linearized upstream solution is then 

nux/p 
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oo 

= I A n 


0 cos nny 

p 


_ 0 1_ 
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sin nuy 



(85) 



For isentropic flow there are two downstream-running waves propagating 
information to the upstream boundary from outside the computational domain and 
one upstream-running wave propagating information from the numerical solution. 
Equation (84) provides one of the lacking pieces of information from outside 
the computational domain; the remaining information is provided by combining 
equations (22), (78) and (79) to give 



1 - Moo 

B " = 1 + C " 

Using the expansions (20) and the Fourier representation (79), the 
coefficients C are determined from 



( 86 ) 



Rl = l C n cos niry = R num - R„ (87) 

1 



where R is the boundary distribution of R obtained from the nonlinear 
num J 

numerical solution and R is given by equation (75). The coefficients A m and 
B n are then obtained from equations (84) and (86). Using the Fourier 
representations (77) and (78), the distributions of 6 and Q on the boundary 
(i.e., the boundary conditions) are calculated according to 



00 

0b = l A n sin nny 

1 

( 88 ) 

00 

Qb = Qoo + £ B n cos nny 

1 
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For the region downstream of the computational boundary (i.e., x>0), the 
exponentially growing component of the solution (35) can be suppressed by 
forcing Kj defined by equation (82) to be zero. This requires the Fourier 
coefficients to be related by 



Cn = B n + A„ 

P 



(89) 



The linearized downstream solution is then 



'Tl‘ 


oo 

= E A n 


~2q 00 M 00 

’ cos nuy 
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_ e i_ 
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sin ntry 



-nnx/p 



(90) 



For isentropic flow there are two downstream-running waves propagating 
information to the downstream boundary from the numerical solution and one 
upstream-running wave propagating information from outside the computational 
domain. Equation (89) provides the lacking information. Using the expansion 
(20) and the Fourier representations (77) and (78), the coefficients A n and B 
are determined from 



00 

6l = £ A n sin nuy = 6 num 

1 

(91) 

00 

Ql = £ B n cos niry = Q num - 

1 



where 0 and Q are the boundary distributions of 0 and Q obtained from 
the nonlinear numerical solution and is given by equations (75). The 
coefficients C n are then obtained from equation (89). Using the Fourier 
representation (79), the distribution of R on the boundary (i.e., the bounda 
condition) is calculated according to 
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Rb = R co + i C n cos nuy 
1 



4 2q 00 M 00 

= Qnum " a « + p 



2 A n cos nuy 
1 
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The second form of this relation requires only the calculation of the 
coefficients A n> 

Many numerical solution algorithms for the Euler equations use so-called 
characteristic far field boundary conditions in which 0 and Q are specified 

oo oo 

at inflow boundaries and R^ is imposed at outflow boundaries. The boundary 
conditions (88) and (92) therefore represent a logical first-order extension 
of the widely-used characteristic (or zero-order) boundary conditions. 

Non-isentropic Conditions - For the case of non-isentropic flow crossing a 
downstream computational boundary, the Fourier series representation (77) for 
0j at the boundary is still valid because 0 is zero at the duct walls. The 
variable o which characterizes the entropy distribution at infinity can be 
represented by 



o = o w + l Dk sin(k-l/2)ny (93) 

1 

where a is the wall value at y = 0. This Fourier series representation 
w 

assumes an even extension of o for 1 < y < 2. The source term component 
defined by equation (55) can then be written as 



2 ^ 

<#> n = —^~r { I (k-1/2) D|< [cos(n+k-l/2)uy + cos(n-k+l/2)ny] 

Y_1 ^ k=l 

(94) 

+ n m£ l Dk [cos(n+k-l/2)ny - cos(n-k+l/2)uy] - 2o w n sin nny} 
k=l 



The procedure for relating o to the entropy distribution at the computational 
boundary and determining the coefficients will be presented later. 



The approximate solution for linearized non-isentropic flow is given by 

equation (61). Since it is only an approximate solution of the approximate 

system (44) and (45), the arguments used in the isentropic analysis for 

choosing K, and K„ to eliminate exponentially growing solution components 
In Zn 

cannot be used. That is, the RHS of equation (51) has been treated as a known 
source term by using a lower-order solution approximation. As a result, the 
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behavior of the source term for large x directly influences the character of 
the solution obtained. However, closer examination of the isentropic analysis 
shows that the same results would have been obtained had the exponentially 
growing components simply been ignored. The remaining constant (i.e., Kj n or 
K^) associated with the decaying solution component could have been 
determined by matching the solution at x=0 with the Fourier series 
representations on the boundary. 



Using this approach and selecting only the decaying component of the 
solution (61), the solution for the downstream region (i.e., x>0) is 
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2qooMoo 
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where 



fn = 



l y 

P 6A n A n y sin nny - — ( 4> n cos nn(y-n) dn 

J o 



hn = 



1 y 

P 6A n A n y cos nny + — ( 4> n sin nn(y-n) dn 

} o 



(95) 



(96) 

(97) 



= nirp 



i 1 sin n " n 



dq 



(98) 



These last three expressions are obtained from equations (59), (62), and (63) 
If the upstream region (i.e., x<0) were to be analyzed for non-isentropic 
conditions, then the other solution component would need to be selected. 
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The solution (95) can be matched to conditions at the computational 
boundary x=0 by letting 



K 



n 



= A + 6A 
n n 



and using the Fourier series representation (77) to give 



( 99 ) 
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6A m = E 7777 ( sin m "y f 't’n sin nii(n-y) dn dy - 
n=l niT J o Jo 



00 r 1 

2p Z A n 6A n F sin muy cos nuy y dy 
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( 100 ) 



Evaluation of the integrals in this expression along with those appearing in 
equations ( 96 ) - ( 98 ) is given in Appendix A. 



All parameters appearing in the solution (95) can be evaluated in terms 

of the Fourier coefficients A and D, . The A coefficients can be determined 

n k n 

in the same manner as for isentropic flow using equation (91), namely, 



00 

Z A n sin nuy = 0 num (101) 

1 

where 0 is the distribution of 0 along the boundary obtained from the 
num 

nonlinear numerical solution. The coefficients can be determined by 
solving the entropy equation (43). Neglecting higher-order terms, its solu- 
tion is 



Si = . 
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where, from equation (93), 
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(k-1/2) D k cos(k-l/2)ny 



( 102 ) 



(103) 
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Combining these two expressions gives 



S-. = I I A n [sin(n+k-l/2)ny + 

Y_1 k=l n=l n 

-(r 1 + <5A n ) x 

si n(n-k+l/2)ny] e 0 



(104) 



Since entropy is convected downstream along streamlines, the distribution 
at the computational boundary, denoted by S num , is known from the numerical 
solution. From the expansions (14) and the definition (37) it follows that 



S = S + S.(0,y) 
num 00 1 

= -2-r (1 - o) + s.(0,y) 

Y-l 1 



(105) 



Introducing the relationship (93) for o and S^O.y) from the solution (104) 
gives 



00 V - 1 

Z Ok si n ( k- 1/2 ) Try = 1 - o w - S num 
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n 00 
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(106) 



l k ~ Hi. a D. [sin(n+k-l/2)-ny + sin(n-k+l/2)ny] 
n n k 

k=l n=l 



The coefficients D k must therefore satisfy 



D m = 2 J 1 [1 - o w - ^ S num ] sin(m-l/2)ny dy 



2 ThT ° k A m-k - 2 1 ° k Am+k ' 1 + 2 1 ^ ° k Ak_ 

k=l k=l k =m+l 



k-1/2 



(107) 



m 



This equation can be solved for each D m by iteration, started by using the 
first integral term as an initial guess. Note that the wall value is known 
from the numerical solution since the wall is a streamline. 
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Combining the definition (7) and the expansions (14) gives the 
first-order relation 



R = Q - T - T 1 
°° 1 



(108) 



Applying this approximation along with the solution (95) at the computational 
boundary x=0, the distribution of R on the boundary (i.e., the boundary 
condition) can be calculated according to 



The approximation (76) for has been used and Q num denotes the boundary 
distribution of Q obtained from the nonlinear numerical solution. This 
expression reduces to the boundary condition (92) if the flow becomes 
isentropic. 



Both upstream and downstream cascade boundary conditions will be derived 
for isentropic conditions. The analysis parallels that for isentropic duct 
flow except that the parameters appearing in the general linearized solution 
(71) are complex quantities. For isentropic conditions, the cascade turns the 
flow without losses. Therefore, all upstream and downstream far field 
quantities (e.g., Q^, R , etc.) can be determined from the downstream pressure 
at infinity, p w . 



% = Qnum " a « + 



2q„Moo 
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(109) 




Cascade Flow 
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At a computation boundary (assumed located at x=0) the perturbation flow 
variables can be represented by the general Fourier series 
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r + E A n cos mry + I A n sin nny 
L 1 1 



Ql = Z B n e 



mny 
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B, 



y + E B n cos nny + E B n sin nny 
* 1 1 
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E C n e 
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+ E C n cos nny + E C n sin nny 
c 1 1 



Ti = E E n e 



i nny 



(113 



= y + E E n cos nny + I E n sin nny 



1 



1 



The boundary conditions will be developed from relationships between the 
Fourier coefficients A n> B n and C n « Fourier analysis provides the additiona' 
relationships (for n = 0, 1, 2, . . .) 



I = Z + Z 
n n -n 



z n * 5 < z „ - 1 V 
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Z n = ' < Z „ - Z -n> 



Z -n ■ I < ! n * 1 Z „> 
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1 



where Z represents any of the coefficients A, B, C or E. Since perturbations 
must vanish at infinity, A q , B q , C o and E q must each be zero. 

Applying the solution (71) at x=0 and using the series expansions (110) 
and (113) it follows that 



A = K. + K, 
n In 2n 



2q 00 M 00 

E n " p ^ (*^2n - ^i n ) 



(115) 



Therefore , 

K ln s + 2qffi„ (1 16) 

k 2" * \ [ft " • Sl E " ] (n7) 

In the region upstream of the computational boundary (i.e., x<0) the 

inflow angle 0^ is known so that t defined by equation (16) is known. To 

suppress the exponentially growing component of the solution (71), the 

coefficients K, and K~ must satisfy 
in Zn 



K ln =0 (n > 0) 



K 2n =0 (n < 0) 



(118) 



This requires the Fourier coefficients to be related by 
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An = En = 2qjC (Bn " Cn) 

The linearized upstream solution is then 
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(119) 



Tl 


-1 


-2q 00 M 00 i 


inn(Aix+y) 00 


2qooMooi 




= £ A n 


P 


e + l A n 


P 


r— 1 ’ 

CD 


— oo 


1 


1 


1 



inu(A2X+y) 
e c (120) 



29 



The parameters X ^ and X^ are defined by equations (72). The component 
corresponding to n=0 has been neglected. This will be discussed in more 
detail below. 



Equation (86) used for the duct flow analysis is again applicable. 

Equation (87) can be generalized to determine the coefficients C n> That is, 

C and C are determined from 
n n 

00 _ oo _ 

£ C n cos nuy + I C n sin nuy = R num - R„ (121) 

1 1 

where R is the boundary distribution of R obtained from the nonlinear 
num 

numerical solution. The coefficient C has been neglected. The coefficients 

o 

A n and B n are then obtained from equations (86) and (119) and the general 
Fourier relationships (114). Using the Fourier representations (110) and 
(111), the distributions of 6 and Q on the boundary (i.e., the boundary 
conditions) are calculated according to 



00 _ oo _ 

6b = 0a, + £ A n cos miy + I A n sin nuy 

1 1 



Qb = Qoo + £ B n cos nuy + £ B n sin nuy 
1 1 



( 122 ) 



Both A o and B q have been set to zero in the boundary conditions (122). A zero 

value of § Q imposed on the numerical solution through the boundary condition 

should tend to force a zero value for C (and hence E ) in the numerical 

0 0 

solution according to equation (86). It should be noted from these 
first-order relations that, within the context of a Fourier representation, Q m 
and 6 can be viewed as the mean values of Q and 6 on the computational 

OO 

boundary. Deviations from this condition can be shown to be a second-order 
effect. These arguments apply also to duct flows. 
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For the region downstream of the computational boundary (i.e,, x>0), the 
outflow angle 9^ is unknown so that t is also unknown. To suppress the 
exponentially growing component of the solution (71), the coefficients Kj and 
K^ n must satisfy 

K ln = 0 (n < 0) K 2n = ° (n > 0) (123) 

This requires the Fourier coefficients to be related by 



An ~ ’ En = " 2qfe (Bn " Cn) (n 0) 



The linearized downstream solution is then 
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(124) 



(125) 



As for the duct flow analysis, the coefficients A and B are determined 

n n 

from the general Fourier representations (110) and (111) according to 



00 _ oo _ 

Z A n cos mry + z A n sin mry = 0 num - 

1 1 

(126) 

00 _ 00 

Z B n cos mry + Z B n sin niry = Q num - Q 

1 1 

where 9 and Q are the boundary distributions of 0 and Q obtained from 
num x num 

the nonlinear numerical solution. The coefficients S and § have been 

0 0 

neglected. The outflow angle 9^ can be calculated from 
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( 127 ) 



= ( e num dy 
o 

which is the mean value of 0 on the boundary. This justifies setting the 

coefficient A„ to zero. The coefficients C are obtained from equations 
o n ^ 

(124). Using the Fourier representation (112), the distribution of R on the 
boundary (i.e., the boundary condition) is calculated from 
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The second form of this relation requires only the evalution of the 
coefficients A n and A n> It is analogous to equation (92) for duct flow. 



The coefficient C q has been set to zero which, when imposed on the 

numerical solution through the boundary condition, should tend to force a zero 

value for B (and hence E ). That is, the boundary condition tends to enforce 
0 0 

R as the mean value of R on the boundary. This argument applies also to duct 

00 

flows. 
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V. APPLICATIONS 



Numerical solutions of the Euler equations have been calculated for 
two-dimensional, steady duct flows using the first-order boundary condition 
procedures developed in the previous section. Both isentropic and 
non-isentropic cases were analyzed. Non-isentropic flow was produced by a 
shock wave located in a nozzle portion of the duct. First-order boundary 
condition results are compared with those produced using the conventional 
zero-order characteristic boundary conditions. 

A solution algorithm was used for the nonlinear Euler equations (l)-(4) 
which is based on the method presented in Reference 1. It uses explicit time 
integration to relax to steady state conditions and includes the shock-fitting 
procedure of Reference 4 to accurately calculate flows containing shock waves. 
Shock waves appear in the solution as distinct discontinuities which satisfy 
the Rankine-Hugoniot relations. It should be noted that the boundary 
condition analysis is independent of the choice of inviscid, nonlinear 
solution method. 

The duct/nozzle geometry is shown schematically in Figure 1. The flow is 
characterized by p ro , the downstream pressure at infinity, which produces a 
mass flow per unit area w through the duct. The linearized solutions given by 
equations (85), (90) and (95) are assumed valid in the semi-infinite regions I 
and 111 and the computational boundary conditions are applied at the upstream 
and downstream boundaries AA and BB of the nonlinear computational region 11. 

The actual shape of the duct/nozzle and the computational grid are shown 
in Figure 3. The nozzle contour is sinusoidal and symmetric about the 
centerline. The computational grid for this portion of the nozzle had 
dimensions 41 x 21. The area ratio is .75 and the upstream and downstream 
areas are equal. For these constant area sections of the duct, additional 
rectangular grid cells could be added without altering the basic 41 x 21 grid. 
This served to minimize the effect of grid changes on the calculations when 
the computational boundaries were moved in order to assess the accuracy of the 
boundary conditions. 
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Because the configuration is symmetric, calculations were limited to the 
lower half of the duct and a centerline symmetry condition was used. Although 
the configuration used for these calculations is simple, the boundary 
condition analysis of the previous section is general and applicable to 
unsymmetric configurations having unequal upstream and downstream areas. Use 
of the simple configuration is sufficient to demonstrate the validity of the 
analysis. 



Isentropic Conditions 

For isentropic flow the first-order upstream and downstream boundary 
conditions are given by equations (88) and (92), respectively. The associated 
analytic far field solutions are given by equations (85) and (90). The 
zero-order (or characteristic) boundary conditions consist of imposing the 
constant value of and a zero value of 0 along the upstream boundary and the 
constant value of R along the downstream boundary. The values of Q and R 

oo oo oo 

are determined from equations (74) and (75). 

Results are presented for a single value of p^ but with the computational 
boundaries located at several different longitudinal stations. The relative 
accuracy of the zero and first-order boundary conditions can then be 
eval uated . 

Case 1 - Results obtained using the complete grid shown in Figure 3 are 
presented in Figures 4 and 5. This grid has 20 columns of grid cells in both 
the upstream and downstream constant area portions of the duct. For this case 
the computational boundaries were far enough upstream and downstream that the 
zero-order and first-order boundary condition results were nearly identical. 
Figure 4 shows pressure and Mach number distributions along the centerline and 
lower wall of the duct/nozzle. Pressure, Mach number, and flow angle contours 
are presented in Figure 5. These results serve as a reference for evaluating 
the accuracy of solutions where the computational boundaries are moved closer 
to the nozzle portion of the duct. 

Case 2 - Results for a shortened duct obtained using the zero-order boundary 
conditions are presented in Figure 6 and 7. There were 5 columns of grid 
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cells in both the upstream and downstream constant area portions of the duct 
for this case. Pressure and Mach number distributions are shown in Figure 6 
and contours are shown in Figure 7. The most noticeable effect of the 
boundary proximity is a change in peak values of pressure and Mach number and 
a significant amount of longitudinal asymmetry. 

Case 3 - The previous shortened duct case was recalculated using the 
first-order boundary conditions. These results are presented in Figures 8 and 
9. Pressure and Mach number distributions are shown in Figure 8 and contours 
are shown in Figure 9. The results within the numerical solution portion of 
the domain are almost identical to those shown in Figures 4 and 5. Linearized 
solution results obtained from equations (85) and (90) have been added 
upstream and downstream of the computational boundaries. It is evident that 
the linearized far field analytic solutions provide for a smooth transition 
across the computational boundary to the true far field conditions at 
infinity. 

Case 4 - Results for a more drastic case using the first-order boundary 
conditions are presented in Figures 10 and 11. Only one column of grid cells 
was used in the upstream and downstream constant-area portions of the duct. 
Pressure and Mach number distributions are shown in Figure 10 and contours in 
Figure 11. Even with the computational boundaries extremely close, the 
numerical solution in the nozzle portion of the duct is affected very little. 

Boundary Comparison - A more quantitative comparison between the zero- and 
first-order boundary conditions can be obtained by examining the distribution 
of flow variables along the computational boundaries. Figure 12 shows a 
comparison of the distribution of pressure and flow angle along the upstream 
and downstream boundaries for Cases 2 and 3. Results taken from the numerical 
solution of Case 1 at the same longitudinal locations are also shown. There 
is little difference between the Case 1 and Case 3 results. The larger 
differences between the Case 2 and Case 3 pressures occurring along the 
upstream boundary are primarily due to the fact that a zero flow angle is 
imposed along the boundary. 
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Implementation of the first-order boundary conditions obviously 
introduces additional computational effort. However, the boundary conditions 
do not have to be updated with every iteration of the numerical solution. The 
small additional effort is more than offset by the large reduction in number 
of grid points required and the fewer solution iterations necessary for 
convergence because information propagates between the computational 
boundaries more quickly. For the isentropic results presented above, only 5 
Fourier modes were required to give a reasonably accurate representation of 
the boundary distributions of R (upstream) and 0 (downstream) for the Case 3 
calculations. For the Case 1 and Case 4 calculations, the number of Fourier 
modes required were 3 and 7, respectively. 

Non-Isentropic Conditions 



Non-i sentropic flow was produced by lowering the value of p^ so that a 
shock wave formed in the nozzle portion of the duct. For this situation the 
isentropic, first-order upstream boundary conditions (88) and the associated 
linearized solution can be used. The upstream far field quantities Q and R 

oo oo 

are constant and must be determined from the mass flow through the nozzle as 
explained when deriving equations (75). The mass flow can be determined by 
numerically integrating the mass flux crossing the downstream computational 
boundary each time the boundary conditions are updated. 

The first-order downstream boundary condition is given by equation (109) 
and the associated analytic solution by equation (95). It should be noted 
that this solution is only an approximation to the solution of the linearized 
equations ( 23 ) - ( 25 ) . The downstream far field quantities and R^ vary in 
the streamline-normal direction and can be determined from equations (76). 
Zero-order downstream boundary conditions as used herein consist of imposing a 
distribution of R^ along the boundary given by equations (76) with o replaced 
by 1- S num . This in itself represents an improvement over some 
applications of zero-order characteristic boundary conditions in which a 
constant value of the Riemann variable associated with upstream-running waves 
is imposed along the boundary. 

Results are presented for two different values of p^ chosen such that a 
weak and a strong shock are produced in the nozzle. The strong shock spans 
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the width of the nozzle while the weak shock does not. Two computational 
boundary locations are used for each value of so that the accuracy of the 
non-isentropic boundary conditions can be evaluated. 

Case 5 - Weak shock results obtained using the complete grid of Figure 3 are 
shown in Figures 13 and 14. Pressure and Mach number distributions are shown 
in Figure 13 for the centerline and lower wall. Pressure, Mach number, 
entropy and flow angle contours are presented in Figure 14. The non-constant 
Mach number in the downstream far field is evident in Figure 13. Note that 
the scales used in Figure 13 are much more compressed than those used for the 
isentropic results. 

For this case, there were slight differences produced by the choice of 
boundary conditions, so these two sets of results are shown separately. These 
differences can be attributed to difficulties with the numerical fitting of a 
weak shock discontinuity. The shock fitting procedure is extremely sensitive 
to small disturbances in the numerical solution when the shock is weak, 
especially in the vicinity of the tail of the shock located near the center of 
the computational domain. These interactions were stronger when the 
zero-order boundary conditions were used. 

Placement of the computational boundary even further downstream by using 
additional columns of grid cells caused very little change to the solution 
when the first-order boundary conditions were used, but produced noticeable 
differences when the zero-order conditions were used. For this reason the 
results of Figures 13b and 14b can be considered more accurate and can serve 
as a reference for evaluating the accuracy of solutions where the 
computational boundaries are closer to the nozzle portion of the duct. 

Case 6 - Results for a shortened duct using the zero-order boundary conditions 
are presented in Figures 15 and 16. There were 5 columns of grid cells in 
both the upstream and downstream constant area portions of the duct for this 
case. Pressure and Mach number distributions are shown in Figure 15 and 
contours are shown in Figure 16. The effect of the boundary proximity is 
clearly evident in these results. 
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Case 7 - The short duct flowfield was recalculated using the first-order 
boundary conditions. These results are presented in Figures 17 and 18. 
Pressure and Mach number distributions are shown in Figure 17 and contours in 
Figure 18. The results agree closely with those shown in Figures 13b and 14b 
for the numerical solution portion of the domain. Linearized solution results 
obtained from equations (85) and (95) have been added upstream and downstream 
of the computational boundaries. Even though the solution (95) is approxi- 
mate, having been obtained by linearizing the thermodynamic relation (36) and 
using only a two iteration approximation for the solution of equation (49), it 
still leads to a useful boundary condition and an adequate prediction of the 
downstream flowfield. 

Case 8 - Strong shock results obtained using the complete grid of Figure 3 are 
presented in Figures 19 and 20. Pressure and Mach number distributions are 
shown in Figure 19 and contours in Figure 20. The significant variation of 
Mach number in the downstream portion of the duct is clearly evident in Figure 
19. The shock fitting procedure was much less sensitive to numerical 
disturbances for this strong shock case. As a result, the solutions obtained 
using the zero-order and first-order boundary conditions for the long duct 
were nearly identical. These results can serve as a reference for evaluating 
the accuracy of solutions where the computational boundaries are closer to the 
nozzle portion of the duct. 

Case 9 - Results for the shortened duct using the zero-order boundary 
conditions are presented in Figures 21 and 22. Pressure and Mach number 
distributions are shown in Figure 21 and contours in Figure 22. The 
inaccuracy of the zero-order boundary conditions for this case is evident. 

One noticeable effect is the change in shock position. 

Case 10 - The short duct flowfield was recalculated using the first-order 
boundary conditions. These results are presented in Figures 23 and 24. 
Pressure and Mach number distributions are shown in Figure 23 and contours in 
Figure 24. For the numerical solution portion of the domain, the results 
agree closely with those shown in Figures 19 and 20. Linearized solution 
results obtained from Equations (85) and (95) have been added upstream and 
downstream of the computational boundaries. The transition across the 
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boundary is reasonably smooth indicating that the approximations used in 
developing the non-isentropi c boundary conditions do not produce a strong 
mismatch between the numerical and analytic solutions. 

Boundary Comparison - A more quantitative comparison between the zero- and 
first-order downstream boundary conditions can be obtained from the 
distribution of flow variables along the boundary. Figure 25 shows a 
comparison of the pressure and flow angle distributions along the downstream 
boundary for Cases 6 and 7 (weak shock) and Cases 9 and 10 (strong shock). 
Results taken from the numerical solutions of Case 5 (first-order boundary 
conditions) and Case 8 at the same longitudinal location are also shown. The 
validity of the first-order boundary condition analysis is clearly evident. 

The minor deviations are due in part to the approximations used in obtaining 
the solution (95). 

For the non-isentropic calculations, 9 Fourier modes were required to 
give a reasonably accurate representation of the downstream boundary 
distributions of 9 and o for the Case 7 and 10 calculations. For the Case 5 
and 8 calculations, only 7 modes were used. Implementation of the first-order 
boundary conditions for non-isentropic situations clearly requires more 
computational effort than for isentropic flow. However, this increase is 
warranted in view of the fact that the far field flow can contain significant 
streamline-normal variations in quantities other than pressure (e.g., see 
Figure 19). 
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VI. SUMMARY 



Far field computational boundary conditions have been developed for 2D 
duct and cascade flows. These first-order boundary conditions are derived 
from analytic solutions of the linearized Euler equations and represent a 
logical extension of the zero-order (or characteristic) boundary conditions 
commonly used in the numerical solution of nonlinear fluid dynamic equations. 
The boundary conditions and analytic solutions provide a smooth transition 
across a computational boundary to the true far field conditions at infinity. 
The boundary procedure is general in that it can be used in conjunction with 
any numerical solution method. 

For the case of isentropic flow the linearized Euler equations are solved 
exactly for duct and cascade flow using separation of variables and Fourier 
analysis. For non-isentropic flow the far field can contain significant 
streamline-normal variations in quantities other than pressure and flow angle, 
and these variations cannot be treated as small perturbations from uniform 
flow conditions. For the case of duct flow the non-isentropic linearized 
equations are solved approximately using an approximate thermodynamic relation 
and an iterative technique. The validity of this approximate solution has 
been verified by numerical results. Use of non-constant zero-order boundary 
conditions, which vary on the boundary as a result of entropy variations, also 
offer substantial improvement over constant conditions derived solely from the 
far field pressure. Extension of the non-isentropic analysis to cascade flows 
is straightforward. 

Use of zero-order (or characteristic) boundary conditions requires that 
the computational boundaries be located far from the nonlinear region of the 
flow. Closer placement of the boundaries may result in a significant amount 
of solution degradation. The first-order boundary conditions allow the 
boundaries to be located much closer thereby reducing the number of grid 
points needed for the numerical solution and also the number of iterations for 
solution convergence. This allows a significant reduction in the amount of 
computational effort required for the nonlinear numerical solution because the 
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additional calculations required for the first-order boundary conditions is 
modest. 

The boundary condition procedures developed herein can be extended to 
axisymmetric, three-dimensional, and viscous flows, including external flows. 
It is recommended that these analyses be undertaken and also that the present 
results be extended to non-isentropi c cascade flows. 



Acknowledgement - This work was partially supported by Naval Air System 
Command (G. Derderian - Program Manager) under contracts N62271-84-M-1857 and 
N62271-86-M-0202, from the Naval Postgraduate School, Monterey, Ca., where the 
technical monitor was Professor Raymond P. Shreeve. 

The author would also like to thank Mr. David S. Stookesberry for his 
help in performing the computations. 



41 



VII. REFERENCES 



1. Verhoff, A. and O'Neil, P. J., "A Natural Formulation for Numerical 
Solution of the Euler Equations," AIAA Paper No. 84-0163 Jan. 1984. 

2. Verhoff, A., "Modeling of Computational and Solid Surface Boundary 
Conditions for Fluid Dynamics Calculations," AIAA Paper No. 85-1496 CP, 
July 1985. 

3. Hirsch, Ch., and Verhoff, A., "Far Field Numerical Boundary Conditions 
for Internal and Cascade Flow Computations," AIAA Paper, To Be Published, 
1989. 

4. Moretti , G., "A Technique for Integrating Two-Dimensional Euler 
Equations," Computers and Fluids, Vol . 15, No. 1, March 1987, pg. 59-75. 



42 



III. APPENDIX A: EVALUATION OF INTEGRALS 



Evaluation of the integrals appearing in equations (96), (97), (98) and 
(100) will be presented here. They are defined as shown below. 

Equation (96): 

i y 

11 = — ( <J>n cos My-n) d n 

J o 

Equation (97): 

i y 

12 = rvn ( < *’ n S ^ n dr > 

J o 

Equation (98): 

‘3 * {‘ ♦» si " "** »» 

Equation (100): 

I4 = f 1 sin muy cos niry y dy 
■^o 

1 y 

15= f sin muy f <J> n sin nu(n~y) dn dy 

^ 0 J 0 

The quantity $ n appearing in the integrands is defined by equation (94). 
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Evaluation of Ij 
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Evaluation of I 2 
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Evaluation of I3 
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Figure 1. Duct/Nozzle Schematic 






Figure 2. Cascade Schematic 




Figure 3. Duct/Nozzle Geometry and Grid 
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Figure 5. Pressure, Mach Number and Flow Angle Contours 
Zero- and First-Order Boundary Conditions 
p = 0.90 
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Figure 6. Pressure and Mach Number Distributions 
Zero-Order Boundary Conditions 
p = 0.90 
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Figure 7. Pressure, Mach Number and Flow Angle Contours 

Zero-Order Boundary Conditions 
p = 0.90 
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Figure 8. Pressure and Mach Number Distributions 
First-Order Boundary Conditions 
p =0.90 
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Figure 9. Pressure, Mach Number and Flow Angle Contours 

First-Order Boundary Conditions 
p = 0.90 
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Figure 10. Pressure and Mach Number Distributions 
First-Order Boundary Conditions 
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Figure 11. Pressure, Mach Number and Flow Angle Contours 

First-Order Boundary Conditions 
p = 0.90 
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Figure 12. Pressure and Flow Angle Distributions 
on Computational Boundaries 
p = 0.90 
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Figure 13a. Pressure and Mach Number Distributions 
Zero-Order Boundary Conditions 
p = 0.83 
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Figure 13b. Pressure and Mach Number Distributions 
First-Order Boundary Conditions 
p = 0.83 
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Figure 14a. Pressure, Mach Number, Entropy and 
Flow Angle Contours 
Zero-Order Boundary Conditions 

p = 0.83 
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Figure 14b. Pressure, Mach Number, Entropy and 
Flow Angle Contours 
First-Order Boundary Conditions 
p =0.83 
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Figure 15. Pressure and Mach Number Distributions 
Zero-Order Boundary Conditions 
p = 0.83 
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Figure 16. Pressure, Mach Number, Entropy and 
Flow Angle Contours 
Zero-Order Boundary Conditions 

p = 0.83 
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Figure 17. Pressure and Mach Number Distributions 
First-Order Boundary Conditions 
p = 0.83 
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Figure 18. Pressure, Mach Number, Entropy and 
Flow Angle Contours 
First-Order Boundary Conditions 
p = 0.83 
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Figure 19. Pressure and Mach Number Distributions 
Zero- and First-Order Boundary Conditions 
p =0.78 
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Figure 20. Pressure, Mach Number, Entropy and 
Flow Angle Contours 

Zero- and First-Order Boundary Conditions 

p =0.78 
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Figure 22. Pressure, Mach Number, Entropy and 
Flow Angle Contours 
Zero-Order Boundary Conditions 
p =0.78 
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Figure 23. Pressure and Mach Number Distributions 
First-Order Boundary Conditions 
p =0.78 
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Figure 24. Pressure, Mach Number, Entropy and 
Flow Angle Contours 
First-Order Boundary Conditions 
p =0.78 
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Figure 25. Pressure and Flow Angle Distributions 
on Downstream Computational Boundary 
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